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Abstract 

We describe an infinite dimensional family of flows E with the prop- 
erty that the associated dynamical system 

*(<) = £(*(*)). ’ i(0)eR" 

is explicitly integrable in closed form. These flows E are of the form 
E = E\ + £2, where E\ and E2 are the generators of a nilpotent Lie 
algebra, which is either free, or satisfies some relations “at a point.” 
These flows can then be used to approximate the flows of more general 
types of dynamical systems. 


(NAS A-CR-1 84967) VECTOR FIELDS AND 

NILPOTENT LIE ALGEBRAS (California Urrv«) 

20 p CSCL 12A 

G3/64 


*The author is a National Science Foundation Postdoctoral Research Fellow. 
f The author is supported in part by NASA- Ames grant NAG2-513. 


/A/ - CV'C- O 

JZ/7 ^ 6C> 


N89-26624 

Unclas 

0217700 


1 


1 Introduction 

This paper is concerned with the efficient symbolic computation of flows of 
ordinary differential equations of the form 

x(t) = F(x(t)), x(0) € K n . 

In other words, we want to write the corresponding trajectory t — *■ x(t) 
in closed form. Of course, since this cannot be done for most differential 
equations, it is important to find a means of approximating a general flow 
by a flow which is explicitly integrable in closed form. In this paper we 
describe an infinite dimensional family of explicitly integrable flows F; in a 
companion paper [6] we give an algorithm which, given a general flow F, 
will, in some cases, find a flow E which is explicitly integrable in closed 
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form and is close to the original flow F. More generally, one could imagine 
approximating a flow by a flow which is explicitly integrable in closed form, 
flowing for a time step, finding a new approximating flow, and continuing. 

The type of integrable flows discussed in this paper arise in a very simple 
fashion. We describe the basic idea here and refer to the paper for the details. 
We assume that the differential equation 

*(0 = E(x{t)), *(0) € R*, 


has the form 


x 1 (t) = d\ 

i2(t) = a 2 (xi) 

X 3 (t) = a 3 (x i,x 2 ) 

if v(0 = aiv(xi,x 2 ,...,a:Ar-i) 

for polynomials ai,...,a^/. Since the first equation can be integrated in 
closed form to yield Zj as a function of time, and the j th equation depends 
only upon the first through (j — l)th equations, the entire system can be 
explicitly integrated in closed form. 

In many examples, the vector field E defining the flow has some addi- 
tional structure so that it can be written as the sum of two or more pieces. 
For example, suppose E depends upon two controls t — ► ui(£) and t — ► u 2 (t) 
and can be written as E = u\(t)Ei + u 2 (t)L Alternatively, suppose that 
E can be written as the sum of two pieces E\ -f ei? 2 , where 6 is a parame- 
ter. Observe that if E\ and E 2 are vector fields with polynomial coefficients 
which are homogeneous of degree 1, then the flow generated by E is of the 
form above and, hence, is explicitly integrable in closed form. See §2 and [4] 
for more details. 

Section 2 provides some background on vector fields, flows and Lie al- 
gebras. Section 3 gives an algorithm for producing flows associated with 
free, nilpotent Lie algebras. Section 4 considers flows which are modeled 
by nilpotent Lie algebras satisfying relations “at a point”. The appendix 
contains a brief description of the MAPLE code which we used to construct 
the examples. 

Nilpotent approximations of this type have been used in control theory 
by Crouch [1] and [2], Hermes [8], [9], and [10], and Krener [13]; and in 
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partial differential equations by Folland [3], Hormander [11], Rockland [14], 
and Stein [15]. 

2 Vector fields, Flows and Lie Algebras 

In this section, we define homogeneous derivations and the flows they gen- 
erate. One of the interests of homogeneous derivations is that their flows 
are integrable explicitly in closed form. Homogeneous derivations may be 
classified by the algebraic properties of the (Lie) algebras they generate. In 
the next section, we give an algorithm which yields homogeneous derivations 
which generate a Lie algebra of special type, called a free, nilpotent Lie al- 
gebra. In section four, we describe the analogous algorithm for derivations 
which satisfy relations “at a point.” For background and related material, 
see Goodman [5], Jacobson [12], and Varadarajan [16]. 

Let R denote the space of polynomial functions on R^ and let E \ , . . . , Em 
denote vector fields on R yV with polynomial coefficients; that, is, Ei is of the 
form 

Ei = Z°iDj, 

j = 1 

where a] = a^(xi , . . . , xyv) G R are polynomials and Dj = d/dxj. These are 
also called derivations , since they are derivations of the ring of polynomials 
R. Recall that a derivation E of a ring R is a map E : R — * R satisfying 
the identities 


E{a + b) = E(a) + E(b) 

E(ab) = aE(b) + bE(a ), 

for elements a and b £ R. The flow generated by the derivation E = 
Y^jLi a ^j} initial point x° — (x®, ...,x%) G R^ is defined to be 
the solution of the initial value problem 

= a\ (*^i ? • • • , % at) 

•^Tv(^) “ U/v (^1 , • • • , X ), 

xi(0) = x?, xn ( 0 ) = x ° n . 

Let t — ► x(t) denote the solution to this initial value problem. We say 
that the point x 1 £ is reached in time s in case x 1 = x(s). Notice that 
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flows form a semigroup; that, is, if we can flow from x° € to the point 
x 1 , and from x 1 to the point x 2 , then we can flow from x° to the point x 2 . 

We are interested in homogeneous derivations of weight ra. We define 
these informally; for a more formal definition, see [5]. First, assign weights 
to the coordinate variables: say xi is of weight pi, . . xjv is of weight pjv. 
Second, assign weights to monomials, by defining the weight of the monomial 
to be 

wt(x^ • • • ) = Pi ei + • • • + Pn gn • 

Third, a polynomial is called homogeneous of weight m in case it is a sum 
of monomials, each of weight m. Fourth, define the weight of the coordinate 
derivation to be 

Wt (Dj) = d/dxj = -Pj, 

and the weight of the derivation a(x)Dj to be wt(a(x)) — wt (Dj). Fifth, 
a general derivation E{ = YjL 1 a i^j 1S °f weight m in case it is the sum 
of derivations, each of weight m. Sixth, the weight of the second order 
differential operator is defined as 

wt (a(x)D{Dj) = wt(a(x)) — wt (£),-) — wt (Dj), 

and the weight of higher order differential operators are defined in an anal- 
ogous fashion. 

If E is a derivation of weight 1 on R^, then the flow 

x{t) = E(x(t)), (0) = x° eR N 

is explicitly integrable in closed form. This is easy to prove by induction, 
since the coefficient of d/dxi, for example, depends only upon variables 
which have already been computed. 

Given M derivations E\ , . . . , Em, we define the Lie algebra generated by 
£i, . . . , to be the smallest real subspace of derivations of R which is 
closed under the formation of Lie brackets 

[E,F\ = EF-FE ; 


this is denoted 

g (E x ,..., Em). 


The length of a Lie bracket is defined inductively as follows: the generators 
all have length 1. If E is a Lie bracket of length less than or equal to k and 
F is a Lie bracket of length less than or equal to /, then [E,F] has length 
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less than or equal to k + 1. The algorithms in this paper use the length of 
a Lie bracket to define homogeneous derivations in the following way. A 
vector space decomposition of = V\ ® ■ • • ® V r is defined by setting 

Vi = span {E : E is a Lie bracket of length < /}. 

If Cj , . . . , €jf is a basis of VJ, then each coordinate x 3 in the dual basis of co- 
ordinates Xj, . . . , Xjt is defined to have weight /, and homgenous derivations 
are defined as above. 

A Lie algebra is called nilpotent of step r in case any Lie bracket of length 
greater than r is zero. Note that if ,£m are homogeneous of weight 

1, then the Lie algebra they generate is nilpotent of step r and any bracket 
of length / is homogeneous of weight /. This is because if F\ is a derivation 
homogeneous of weight k and F 2 is a derivation homogeneous of weight /, 
then [F \ , F 2 \ is homogeneous of weight k + /. A Lie algebra is free in case it 
satisfies as few relations as possible. Note that a Lie algebra always satisfies 
some relations, since for example 

[Er,E 2 } = -[E 2 ,E 1 ] 

and 

[[E u E 2 ] j E 3 ] + [[&,&], E x ] + [[E 39 E 1 ],E 2 ] = 0. 

Also, there will be additional relations that are a consequence of these. A 
Lie algebra is free in case any relation is a consequence of relations that are 
of the form above. Rather than try to make this precise, we give a basis for 
a free, nilpotent Lie algebra that is due to M. Hall. 

Definition 2.1 Each element of the Hall basis is a monomial in the gen- 
erators and is defined recursively as follows. The generators 
are elements of the basis and of length 1. If we have defined basis ele- 
ments of lengths 1 - 1 7 they are simply ordered so that E < F if 

length(E) < length(F). Also if length(E ) = s and length(F) — t and 
r = s + t, then [E, F ] is a basis element of length r if: 

1. E and F are basis elements and E > F, and 

2 . ifE = [G,H], then F > H. 

3 Free, Nilpotent Lie Algebras 

In this section we provide motivation by reviewing an algorithm, whose 
input consists of a rank r, and whose output consists of two vector fields 
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Ei and E 2 on .which generate a Lie algebra which is isomorphic to the 

free, nilpotent Lie algebra g 2 , r on two generators of rank r. Here N is the 
dimension of the Lie algebra. The proof can be found in [4]. In the next 
section we consider the analogous problem when the Lie algebra satisfies 
relations. Recall that the flows of these vector fields can all be integrated 
explicitly in closed form. 

Fix the rank r > 1 of the free, nilpotent Lie algebra g 2 , r , and number 
the basis elements for the Lie algebra by the ordering from Definition 2.1, 
i.e., E 3 = [E 2 >Ei\, E a = [E 3 ,Ei\, E$ = [i? 3 ,i£ 2 ], etc. Consider a basis 
element E * as a bracket in the lower order basis elements, [E^ > ^3\ h where 
H > k • If we repeat this process with P tl , we get Ek = [ [Ei 2 , Ej 2 ] , Pj, ] , 
where j 2 < j 1 by the Hall basis conditions. Continuing in this fashion, we 
get 

Ek = [[' " [ [Ei m , E Jm ], Ej m _ 1 ], • • • , Ej 2 ], Ej x ], 

where i m =2, j m = 1, and j n+1 < j n for 1 < n < m — 1. This defines a 
partial ordering of the basis elements. We say that Ek is a direct descendant 
of each Ei } , and we indicate this by writing k y ij. 

Define monomials P 2t k inductively by P 2< k = —Xj A,t'/(deg J P 2 ,i + 1), 
whenever Ek — [ E{,Ej ] is a Hall basis element, and where deg^P is the 
highest power of Xj which divides P. If l >z m, then define to be the 
quotient P 2 ,m/p 2 ,;- Note that P !tt = 1. The coefficients of these monomials 
guarantee that — g§- Pi,i — Pk,l whenever i -< k l. For example, if Ek = 
[[[[[[[E 2 ,E 1 ),E 1 lE X l iE 2 },E 4 }’E 4 },E 7 }^nd if P, = [[[E 2 ,E 1 ],E l ],E l ], then 


P2,k = - 


XiX 2 X4X 7 

3!2! ' 


and P{ t k — 


x 2 x\x 7 

2 ! 


The following theorem summarizes the properties of the generators. 


Theorem 3.1 Fix r > 1 and let N denote the dimension of the free, nilpo- 
tent Lie algebra on 2 generators of rank r. Then the vector fields 


Ex 

E 2 


d 

dxx 


_d_ 

dxi 


have the following properties: 


1. they are homogeneous of weight one with respect to the grading 


K N = Vx © • • • © V r , 
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where V t = the span of Hall basis elements of length i. 

2. the Hall basis elements E * they generate satisfy £*(0) = g|^; 

3 . £/ie flow E\ + £ 2 25 explicitly integrable in closed form; 
the Lie algebra they generate is isomorphic to g2 t r* 

4 Vector Field Models With Relations at the Ori- 
gin 

In the previous section, we gave an algorithm which generates the explicitly 
integrable flows associated with free, nilpotent Lie algebras. Our eventual 
goal is to find as large a class as possible of explicitly integrable flows. In 
this section, we generalize the algorithm of the previous section by giving 
an algorithm which generates the explicitly integrable flows associated with 
nilpotent Lie algebras satisfying “relations at a point.” 

Let £\,...,£n denote the basis of the free, nilpotent Lie algebra on two 
generators £\ and £2 of rank r and dimension TV, and let £ 1 , . . . , Ejy denote 
the vector fields obtained by replacing each occurrence of the generators £\ 
and £2 in £ t by £1 and £ 2 . Recall that a homogeneous relation is a relation 
of the form 

N 

£i = ^ f{ j£j, 

3 = 1 

where r,*j = 0, unless wt (£») = wt(£j), and r t *j are scalars. We are interested 
in quotient algebras g defined by imposing homogeneous relations among 
the generators. Our eventual goal is to find a vector field model for an 
arbitrary quotient algebra g; that is, vector fields £ 1 , . . . , £yv corresponding 
to the Lie elements £\ , . . . £m and satisfying the same bracket relations as the 
quotient algebra. The reason is simple: many questions about the structure 
of the algebra and about the associated dynamical system x(t ) = E\(x(t)) + 
£ 2 (x(i)) are reduced to easily computed questions about the flows of the 
£/s. 

In this paper, we solve just a “local” version of this problem. To de- 
scribe this problem, we fix the basis £\ , . . . , £n in the domain, and the basis 
D \ , . . . , Djv in the range, and use the homogeneous relations to induce a map 
R : — * R^. The problem we solve is to find vector fields E\ and £2 such 

that the associated Lie brackets evaluated at the origin £i(0), . . . , £tv(0) are 
the columns of the matrix R . We prove the theorem assuming the following 
technical hypothesis. 
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Assumption. The. map R is the identity when restricted to the subspace 
spanned by the derivations Si of weight less than or equal to r/2. In other 
words, the components of R satisfy , whenever the weights of Si 

and Sj are both less than r/2. Here 6{j is the Kronecker delta. 

Since the number of basis elements grows very quickly with the weight, 
this requirement affects very few of them. This restriction reduces the ar- 
guments essentially to the proof of the free case; under this condition, the 
images of the free nilpotent Lie algebra vector fields E\ and E 2 under the 
map R generate the desired quotient algebra with relations at the origin. 
Indeed, the same map R takes the basis elements for the free nilpotent Lie 
algebra to a spanning set for the quotient algebra. Of course, in general, 
the operations of bracket and non-isomorphic linear mapping do not com- 
mute. With the above condition, it follows that they at least commute at 
the origin. 

We conclude this introduction with a simple remark. Let B denote 
the indices in the index set {1,...,A} which are of weight less than r/2; 
and let A denote the remaining indices. Also, let R = (f z j) denote the 
matrix defining the homogeneous relations satisfied by the generators £ t -, 
and assume that 


1. fij = Sij , for ij e B 

2. for each fixed i 6 A, r ZJ - = 0, for all j € A and j ^ i. 

It is easy to see that, with these hypotheses, the matrix of relations R is 
equal to the matrix R defined above using the bases S x and Dj. 

Define polynomials P 2 y k inductively by P 2 ,* = -Xj P 2 , t / (deg jP 2 ,i + 1), 
whenever Ek = [Ei,Ej] is a Hall basis element, and where deg^P is the 
highest power of xj which divides P. If l >: m, then define Pi t7n to be the 
quotient P 2 ,mlP 2 J- Note that P tjt = 1. 

Theorem 4.1 Choose R = (r t j) satisfying the assumption above , and set 


Pi 


_ 9 _ 

dxi 


and 


N 


t>2 j-1 



Then the Lie brackets Ej have the properties: 
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1. they are homogeneous of weight 1 with respect to the grading 

R n = Vi © • • ■ 0 V r , 

where V{ = the span of Hall basis elements of length i; 

2. the flows of the system 

x(t) = Ei(x(t)) + E 2 (x(t)), x(0) = x° 

are explicitly integrable in closed form; 

3 . at the origin , the vector fields satisfy 

3 - 1 3 

for all 1 < i < N , and so they satisfy the desired relations. 

Proof of theorem. By construction the vector fields are homogeneous of 
weight 1; therefore their flows are explicitly integrable in closed form; hence, 
we need only prove assertion three. As in the proof of the free case, tins 
is proved using a lemma which is actually stronger than the theorem and 
which illuminates the structure of the vector fields E r We must introduce 
the minimum order of a polynomial which was defined in [4]. 

Definition 4.2 Let A be any non-constant monomial in the variables X \ , 
. . . , x n . Define the minimum order of A by m(A) = min{j : Xj \A}. If A is 
a polynomial with no constant term , then define m(A) to be the maximum 
of the minimum orders of the monomials of A. 


Lemma 4.3 If Ek = [Ei,Ej], then 


Ek = Yl E rm ‘^~ + E 

lyk rn UXm n UXn 


where j < an d Qk,n I s a polynomial with no constant term with 

minimum order < j . 


Assertion three follows easily from the lemma. The only coefficients which 
do not vanish at the origin are from the first term when l = k. This implies 


Ek{ 0) = ELi r mkd § as desired. 


10 



Proof of lemma. -The proof is by induction. Direct computation yields 


£> = E p-E'-i 


»>3 


dxj 


since P 2 ,i = — * 1 ^ 3 , » for alii > 3. 

Consider the parents of By the inductive hypothesis: 


Ei ~ £ Pi ' a ? rfca 9x >> + ? Qi,c dx * ’ 


and 


" dTj + £ Pi4 ^ rde dT e + ^ Q}J ~dT s 


since the weight of < p/2, and so £j(0) = 

Now E{ = [jE 7 p , and so q < j by the basis condition. Therefore, we 

have the relations: 

q < Tn(P it a) < i and m(<2, iC ) < q <j. 

And trivially, 

m {Pj,d)<j and 

Examine the terms in the bracket of E{ and Ej. 


[Ei,Ej] ~ £ £ P ’ ,a r6a Ted ( dx b P] ' d } dx e + £ £ Pi ' a Tba ( dx b Qhf } dxj 

dyj e “* f 

+ V^ Qi ' cred( d^c Pj ' d) dT e + ^^ Qi ' c( dT c Qjj) fa 7 

dyj 1 

- ?<£«*>£ 

E E F ’' dT ’ d, 7iz Q '-'- ] nZ 

a>-» b 


E E dXf p ^) dXb E E Qj>f( d d Xf 

a>-i ■* 
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If A and B are any monomials with m(A) < j then for any m, A^^B 
either vanishes or has minimum order less than j. This, and the fact that 
rij = 0 unless E t and Ej have the same weight, imply that the non-zero 
terms in the second, fifth, and sixth lines have minimum orders < j. If A is 
any monomial satisfying m(A) < j and if m > j, then q§^A either vanishes 
or has minimum order less than j. This implies that the non-zero terms in 
the first and fourth lines have minimum orders < j. The remaining terms 



When a = i, this term vanishes. If m(P l?a ) = j, then either k = a, in which 



or a >- k, and 

~ E rha ^ Pi ' a) jT h = E rbaPk ' a ~di' h 

If m(P t>a ) < j, then -Pi, a is either zero, or it has minimum order less 

than j. If m(Pi A ) > j, then g| -Pi, a = 0. We conclude that 

= EE r mi p u-Q- + 

a>i b ° X 3 ° Xb lyk m 0X ‘ g OX 9 

where m(Qk, g ) < j ■ This proves the lemma. 

A Appendix: Maple Code 

The algorithms described in this paper were the result of computer exper- 
imentation using the symbolic packages MACSYMA and MAPLE. In this 
appendix, we give examples of some of the very simple MAPLE code that 
we wrote. Figure 1 gives the MAPLE code to compute Lie brackets, and 
Figures 2, 3 and 4 give the MAPLE code to compute vector fields which sat- 
isfy relations at a point. This code is illustrated by computing two vector 
fields which satisfy the relation 

Ee( 0) + 2E 7 (0) = E s ( 0). 
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brac:= proc(v.w) 

treturns the Lie bracket of the vector fields v and w 
local temp,t p i, j ,1; 

if not type ( v, 'vector ‘ ) then 
ERROR (‘ Not a vector ') 
elif not type (w , * vector ') then 
ERROR (‘Not a vector ') 
f i; 

if vectdim(v)Ovectdim(w) then 
ERR0R(‘ Different dimensions!') 

fi; 

temp:=array(l . . vectdim(v) ) ; 
t :=array(l . . vectdim(v)) ; 

for j from 1 to vectdim(v) do 
for i from 1 to vectdim(v) do 

t[i] := v [i] *dif f (w [ j] p x [i] ) -w [i] *dif f (v [j] ,x[i]) 
od; 

temp [ j] := sum(t [1] , 1=1 . . vectdim(v) ) ; 
od; 

eval(temp) ; 
end; 


Figure 1: Maple code to compute Lie brackets. 

Finally, in Figure 7, we give the code which will flow along vector fields, 
and illustrate it by flowing along the sum of the two vector fields already 
computed. 

Our example is a quotient algebra of g2,4- This has a basis given by 
Ei, £2 (the generators), £3 = [£2, £1], £4 = [£ 3 ,£ i ], £5 = [£3, £2], E& = 
[£ 4 ,£i], £7 = [£4, £2], and £g = [£5, £4]. We impose the relation that 

£ 7 ( 0 ) + 2 £ 8 ( 0 ) = £ 6 ( 0 ). 

The r matrix, followed by the two vector fields and their brackets is in 
Figures 5 and 6. The result of flowing along the £1 + £2 for time r using 
the procedure in Figure 7 is displayed in Figure 8. 
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HallRelations := proc(rho, maxdim) 

#rho is the maximum weight, maxdim is an upper bound for the dimension 
local left .right ,i,j ,a,b,w,wO,wptr ,k,poly ; 

left := array(l . .maxdim, sparse, C(3)=2] ) ; 
right := array (1 . .maxdim, [(3)=1]) ; 

poly: 2 * array (1 . .maxdim, 1 . .maxdim, sparse, [(3,1)=1]) ; 
w:=array(l . . maxdim, sparse , [(1)*1 , (2)*1 , (3)*2] ) ; 
wptr := array (1. .rho+1 , sparse , [( 1) =1 , (2) =3, (3)«4]) ; 

#wptr(i) is the index of the first element of weight i. 


Figure 2: Maple code to compute vector fields satisfying relations at point: Part 1. 
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1 : - 4: 

print (e .3=[e .2 ,e . 1] ) ; 
for wO from 3 to rho do 

#w0 is the weight of the produced brackets 
for i from round(w0/2) to (wO-1) do 

#i is the weight of the left parent 
for j from wptr[i] to (wptr [i+1] -1) do 

is the index of the left parent 
for k from wptr[wO-i] to (wptr [wO-i+l]-l) do 

#k is the index of the right parent 
if k < right [j] then next fi; 

#These are the Hall Conditions 
if j <= k then break fi; 
left [1] := j; 

#assign the left and right parents 
right [1] := k; 

#to the new bracket, 
print (e.l=[e.j ,e.k] ) ; 
w[l] w[j]+w[k] ; 

#This is its weight, 
for m from 1 to k~l do 

#This loop assigns the multiindex for the coefficient, 
poly [l,m] := poly [j ,m] od; 
poly [1 ,k] := poly[j ,k] + l ; 

1 : = eval(l+l) ; 
od; 
od; 
od; 

wptr[wO+l]:= 1; 
od; 


Figure 3: Maple code to compute vector fields satisfying relations at point: Part 2. 
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el:=array(l. .wptr [rho+l]-l , sparse, C (1 ) =1] ) : 
e2 : =array ( 1 . . wptr [rho+l]-l , sparse , C(2)=lj) ; 
for i from 2 to rho do 

#Sum over the weights 
for j from wptr[ij to wptr[i+l]-l do 

#Sum over elements of the same weight 
for q from wptrfi] to wptr[i+l]-l do 
e2[j] := eval(e2[j])+ 

(product ( (~x [p] ) ~poly[q,p]/ (poly[q,p] !) , p=l . .right [q] )*r [j ,q] ) 

od; 

od; 

od; 

print (el=eval( el)) ; 
print (e2=eval(e2) ) ; 
for i from 3 to wptr [rho+1] -1 do 
a:* left [i] : 
b: =right [i] : 
e.i:= brac(e , a,e.b) ; 
print (e . i=eval(") ) ; 
od; 
end; 


Figure 4: Maple code to compute vector fields satisfying relations at point: Part 3. 
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r : * array ( 1 . . 8 , 1 . . 8 


Cl, 

0, 

0, 

0, 

o, 

0, 

0, 

0] 

[0, 

1, 

0, 

0, 

0, 

0, 

o, 

0] 

[0, 

o, 

1, 

0, 

0, 

0, 

0, 

0] 

[0, 

o, 

o, 

1, 

o, 

o, 

0, 

0] 

[0, 

0, 

o, 

o, 

1, 

0, 

o, 

0] 

[0, 

0, 

0, 

0, 

o, 

0, 

0, 

0] 

[0, 

o, 

0, 

0, 

o, 

1, 

1, 

0] 

[0, 

0, 

0, 

0, 

0, 

2, 

0, 

1] 


HallRelations(4,8) ; 

e3 = [e2, el] 
e4 = [e3, el] 
e5 = [e3, e2] 
e6 = [e4, el] 
e7 = [e4, e2] 
e8 = [e5, e2] 

Figure 5: The relation matrix, the vector fields, and their brackets: part 1. 


17 



el = (array ( sparse, 1 .. 8, [1, 0, 0, 0, 0, 0, 0, 0])) 
e2 = (array ( sparse, 1 . . 8, 

2 3 2 

[0, 1, - x[l], 1/2 x[l] , x[l] x [2] , 0, - 1/6 x[l] - 1/2 x[l] x[2] , 

3 2 

- 1/3 x[l] - 1/2 x[l] x [2] ])) 

e3 = (array ( 1 .. 8, 

2 2 2 

[0, 0, 1, - x[l] , - x[2], 0, 1/2 x[l] + x[l] x [2] , x[l] + 1/2 x[2] ])) 
e4 * (array ( 1 .. 8, 

[0, 0, 0, 1, 0, 0, - x[l] - x[2] , - 2 x[l]])) 
e5 = (array ( 1 . . 8, [0, 0, 0, 0, 1, 0, - x[l] , - x[2]])) 

e6 = (array ( 1 .. 8, [0, 0, 0, 0, 0, 0, 1, 2])) 

e7 = (array ( 1 . . 8, [0, 0, 0, 0, 0, 0, 1, 0])) 

e8 ” (array ( 1 . . 8, [0, 0, 0, 0, 0, 0, 0, 1])) 

Figure 6: The relation matrix, the vector fields, and their brackets: part 2. 
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f low:=proc(v,pO) 

#flows returns the point reached after flowing for time tau 
#along the vector field v from the point pO 

local xt ,r ,xl,i, j ,n; 

n := vectdim(v) ; 
xt : = array (1 . .n) ; 
xl := array (1 . .n) ; 
r := array (1 . .n) ; 
r[l] :=v[l] ; 

for i from 1 to n-1 do 

xt [i] := int(r[i] , t)+pO[i] ; 

r [i+1] : - subs( (x[j]=xt[j]) $ j=l . . i,v[i+l] ) ; 
od; 

xt [n] := int (r [n] , t ) +p0 [n] ; 
for i from 1 to n do 

xl[i] := subs(t=tau,xt [i] ) 
od; 

eval (xl) ; 
end; 


Figure 7: Maple code to flow along a vector field. 


flow (add (el ,e2) , array (1 . .8, sparse)) ; 
array ( 1 . . 8 , 

[tau, tau, - 1/2 tau , 1/6 tau , 1/3 tau , 0, - 1/6 tau , - 5/24 tau ]) 
subs(tau=l , ") ; 

array ( 1 .. 8, [1, 1, -1/2, 1/6, 1/3, 0, -1/6, -5/24]) 


Figure 8: The result of flowing along the vector field el + e2. 
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